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Abstract 

Vectorial and anomaly free chiral U(l) fermion models on a 2d finite lattice are 
considered. It is demonstrated both numerically and analytically that introduc- 
tion of Pauli - Villars type regularization supresses the symmetry breaking effects 
caused by the Wilson term. 
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1 Introduction 



Treating chiral fermions on the lattice still presents serious problem. Due 
to well known " no-go" theorem [|I| there is no straightforward way to remove 
fermion spectrum degeneracy without breaking chiral invariance. Several 
possibilities to overcome this difficulty were proposed (for recent review see 
[0]), but all of them have certain problems. A really successful model has to 
provide reasonable results both in perturbative and nonperturbative regions 
at least for sufficiently small lattice spacing. 

Existing computer facilities allow efficient nonperturbative tests only in 
two-dimensional models, the U(l) gauge model on the torus being a popular 
example. 

In the present paper we apply to the toron model the method, proposed 
in ref.0]. The idea of this method is to introduse in the lattice models 
additional gauge invariant Pauli - Villars (PV) type regularization which 
supresses the contribution of momenta close to the border of Brillouin zone, 
\p\ ~ 7r/a. Pauli - Villars fields mass M introduces a new scale which is 
choosen to satisfy the condition M <C a^^. With this condition fulfilled 
any modification of the action at the distances of order a becomes irrele- 
vant. Therefore, if the additional regularization respects vectorial and chiral 
gauge invariance, possible symmetry breaking effects due to introduction of 
Wilson term or some other device removing spectrum degeneracy, are su- 
pressed and vanish in the continuum limit. 

In the paper it was shown that introducing PV type regularization 
together with the standard Wilson term for anomaly free chiral gauge 
models on the infinite lattice one gets in the framework of perturbation 
theory correct continuum results without chiral noninvariant counterterms. 
Possible chirality breaking effects are of order a. 

Below we shall check these results nonperturbatively for the two dimen- 
sional model - U(l) gauge invariant interaction of fermions on the 2d finite 
lattice. To compare our calculation with known exact results for continuum 
theory on the torus [5 - 7] the gauge field is chosen to be constant. We 



2 



present the results both for vectorial and anomaly free chiral models and 
compare them with the corresponding calculations for Wilson fermions sup- 
plemented by gauge noninvariant counterterms as proposed by the Roma 
group H. 



2 Vectorial lattice model 

We start with the vectorial model on the finite lattice described by the 
action 



hw = \y. [^{x)-i^U^il){x + A) - i^ix) [^{x + ft) - ip{x)]} + h.c. (1) 

Z X,fJ, 

Here -N/2 + l<x^,< N/2, /i=0, 1; is the number of lattice sites. The 
lattice spacing is chosen to be equal to 1. The first term describes gauge 
invariant interaction of fermions with the constant field 

The second term is the Wilson term removing fermion spectrum degeneracy. 
Being interested in the study of symmetry breaking effects on a lattice we 
choose the gauge noninvariant Wilson term. Note that chiral invariance is 
broken both for noninvariant or covariant Wilson terms. The Fermi field ijj 
satisfy antiperiodic boundary conditions. 

A straightforward calculation gives the following expression for the de- 
terminant: 

"^^-.^ J/2.1 5^(P.0) + H/^(P)' 



where 



27r 

B^{p, h) = sin — {pf, -h,- 1/2) , (3) 



I 2tt 
TV 



W{p) = Y. (l-cos — (j)^-l/2) . (4) 



^=0 
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The determinant (|2D is normalized to 1 at h=0. The corresponding expres- 
sion in the continuum theory looks as follows [5-7]: 

Dvc = e-'"'^^ n \F[n, h]F[n, -h]\\ (5) 

n=l 

Here 

^ _^ g-27r(n-l/2)+27ri(/io+i/ii) 
^[^' ^] = 1 ^ g-27r(n-l/2) " 

The determinants (0), (§) satisfy the following symmetry properties: 

D[ho, hi] = D[hi, ho] = D[-ho, hi] = D[ho, -hi]. 

It allows to consider the fields h^ in the interval 0<ho<hi. The continuum 
determinant (g) satisfies also periodicity condition 

Dvc[ho, hi] = Dvc[ho + no, hi + ni], no, rii = 0, ±1, ±2, . . . 

Due to the breaking of gauge invariance by the Wilson term, the lattice 
determinant satisfies weaker periodicity condition 

N N 
Dvw[ho, hi] = Dvw[ho + yno, hi + yni]. 

Computer simulations were performed for0</io<0.5, — 0.5</ii<1.5 us- 
ing the lattices with N=32, 160. The results are presented at Fig.l (curves 
1, 2). One sees that the lattice results do not agree with the continuum 
ones in the whole interval of ho, hi. It shows that although formally the 
action (|I]) has the correct continuum limit, symmetry breaking effects due 
to noninvariant Wilson term change quantum corrections drastically. 

Below we show that this discrepancy vanishes if one modifies the action 
(H]) by introducing additional PV regularization. The modified action looks 
as follows: 

Ipv = ^ E|E \^{x)-f^U^,(l){x + A) - + ft) - + 



-M(J){x)(J){x)\ +h.c 
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Fig.l. Vectorial determinants Dy as functions of hi at ho = 0.2: 

1 - Dvc on the torus; 2 - Dyw without PV field, N = 160; 
3, 4 - DvR with PV field, M = M2: 3 - = 32, 4 - 160; 
5, 6 - DvK with counterterm: 5 - = 32, 6 - 160 



Here is a bosonic PV field having the same spinorial and internal 

structure as tjj. In the considered model one PV field is sufficient to supress 
the contribution of the region near the border of the Brillouin zone. 
Regularized determinant may be presented in the form 

DvR = Dvw[h]Dpv[h], 

where Dvw[h] is given by eq.(@) and Dpy is the corresponding expression 
for the PV field: 



/ B^p, h)^{W{p)^Mf \ ' 



Dpv[h]= n i z.ririii^'rzu ] ■ (6) 
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Fig. 2. Vectorial determinants Dy as functions of M 
at ho = 0.2 and hi = 0.4: 

1 - Dye on the torus; 2, 3 - D^/i? with PV field: 
2- N = 160, 3 - 320 



Note that in this equation we defined the determinant of a bosonic PV 
field as inverse of the corresponding fermion determinant. One can won- 
der if this determinant may be really presented as a path integral of the 
corresponding bosonic action Ipv- A naive integral over bosonic fields of 
eyip{—Ipv) does not exist as the action is not positive definite. However 
as eq.(0) shows the fermionic determinant Dyw is positive and can be 
presented as a path integral of an exponent of a positive definite action. 
Therefore the corresponding PV determinant also can be written as a path 
integral of an exponent of a positive action. Its explisit form can be easily 
read of the eq.(0). In general, if the modulus of a fermion determinant 
can be presented as a path integral of exponent of a positive local action, 
analogous representation is valid for the corresponding PV determinant. 

The results of calculations for regularized determinant as function of the 
regularizing mass M for fixed ho, hi at A^=160 and 320 are presented at 
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Fig. 2. Exact agreement with the continuum is achieved for two values of 
M=Mi, M2. However in the whole interval Mi < M < M2 which in our case 
is of order Mi=0.01 - 0.03 and M2~0.3 the discrepancy is within 10%. So 
to get the correct result one needs to tune the regularizing mass, but it is 
not really a fine tuning. 

The dependence of the regularized determinant on the hi at /io=0.2 for 
value M=M2 are shown at Fig.l (curves 3, 4). One sees that for —0.5 < 
hfj,<0.5 there is a very good agreement with the continuum results both for 
N=32 and 160. To get a good agreement for larger h one needs a bigger 
lattice: for 0.5</ii < 1.5 A^=160 provides a good agreement whereas A^=32 
gives a sizable mistake. 

The reason for that is easily understood. PV regularization supresses the 
contribution of momenta for which B^^W^':^ M^. One sees from the eq.(0) 
that depends only on the difference p^ — h^. Hence for h big enouth B^ 
can be much bigger than even for relatively small p. In this way part 
of a physical region will be also supressed. This effect obviously disappears 
with increasing of N. It seems to be a common problem for any non gauge 
invariant lattice scheme. In particular our calculations show that the same 
phenomenon occurs in Roma approach In the particular model one 
can get better agreement by using other approaches e.g. overlap formalism 
[p] or staggered fermions as was done in somewhat different context in 
ref. 0. In our case we could easily avoid the problem by using a gauge 
invariant Wilson term. In this case our regularized action is manifestly 
gauge invariant and the determinant is periodic with the period 1. Hence 
if there is an agreement with the continuum in the interval —l/2<h^<l/2 
it is continued automatically to any h. 

However for general chiral models considering large nonsmooth external 
fields really makes a problem. One way to avoid it is to supress the con- 
tribution of such fields by some additional regularization. In the paper 
it was shown that in perturbation theory higher covariant derivative for 
gauge field may play this role. It would be important to check if it also in 
nonperturbative calculations. 
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Fig. 3. Values —6v and Q.26c as functions of hi + hi: 
1- N = 160, ho = and 0.2; 2 - 32, 



The alternative way to improve the agreement with the continuum results 
is, following the Roma approach to introduce to the action (|T]) a gauge 
noninvariant counterterm. Introducing the counterterm ky [hi + hfj , one 
gets for the determinant 



To determine the coefficient ky we calculated the value Sv = — ln{Dvw / Dye) 
as a function of hi + hi, the result being presented at Fig. 3. It gives 
/cv" = —1.6877. Substituting this value to eq.([7|) one gets for Dyx the curves 
presented at Fig.l (curves 5, 6). Qualitatively we have the same situation 
as in the case of PV regularization: for N=32 there is a good agreement in 
the interval — 0.5 < /i^ < 0.5, for N = 160 the interval of agreement extends 
to -0.5</i^<1.5. 

Therefore in the model under consideration both PV regularization and 




(7) 
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Roma approach give analogous results. The advantages of PV approach are 
twofold. First of all in this case the only source of gauge symmetry breaking 
of the action is the Wilson term and no gauge noninvariant counterterms 
are needed. One can avoid gauge symmetry breaking completely by choos- 
ing other methods of spectrum degeneracy removing. SLAC-model sup- 
plemented by PV regularization is the example of manifestly gauge invari- 
ant (although nonlocal) model without spectrum doubling [jTD|. Moreover, 
gauge symmetry preserving PV type regularization may be constructed also 
for anomaly free chiral models [0, [TO], [TTj| . 

Secondly, in more complicated four dimensional models in the framework 
of Roma approach one would need to introduce more gauge noninvariant 
counterterms, whereas PV type regularization cuts all unwanted effects of 
momenta p ~ 7r/a, provided the Yang - Mills fields are also regularized 
by means of higher covariant derivatives (see the discussion in [^). It is 
worthwhile to mention that even in a two-dimensional models a simple 
mass renormalization may be not sufficient to get a good agreement with 
the continuum [^. 

The numerical results obtained above are in a good agreement with the 
analitical estimaties. One can show that all the lattice diagrams with more 
than two external lines differ from the corresponding diagrams in the con- 
tinuum toron model by the terms of order 0(1/M7V). At the same time 

-,n 2 

n^^(o) = Uycio) - 

where 11^(7(0) = 27r is the mass gap in the comtinuum toron model. There- 
fore to get the correct continuum result without PV regularization one 
has to add to determinant Dyw the counterterm with the coefficient ky = 
—87r^/27\/3« —1.6884, to be compared with the numerical value —1.6877. 

If the PV fields are introduced the following expressions for polarization 
operator can be derived: 

Uvr{0) = 27r + 0(1 /MN) + 0(M In^ iV), if M ^ when ^ oo; 

nFi?(0) 27r, if M = 0.307 when A^ ^ oo. 
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These values are in a good agreement with the numerical results presented 
at Fig. 2. Analogous estimates are done for anomaly free chiral models. 

3 11112 lattice model 

In this section we apply the approach described above to anomaly free 
chiral models on the finite lattice. The possibility to supress the contribu- 
tion of momenta close to the border of Brillouin zone in anomaly free 4d 
models by introducing a chiral gauge invariant PV type regularization is 
related to the fact that in this case divergent diagrams (with less than 5 
external vector lines) contribute only to the modulus of the determinant. 
From the point of view of these diagrams the theory is essentially vector 
like. In two dimension the anomalous diagram is the polarization operator, 
and the only way to get rid off anomaly is to use a combination of left- 
handed and right-handed interactions in which the terms proportional to 
75 cancel. In particular one can consider 11112 model including four left- 
handed fermions with charge 1 and one right-handed fermion with charge 
2. 

The action of this model can be written in the form 

1 4 . _ _ 

lew = o E E V^fc+(a^)7M {P+U^i + P ) i^k+{x + A) - iJJk+{x) X 

^ k=l ^ x,fi 

x{tjjk+{x + jl) - tlJk+{x))]] + lT.['^-{x)-f,{P+ + Uf,P^) X 

xtj;_{x + /i) + i/j_{x){i/j^{x + /i) - iIj-{x))] + h.c. (8) 

Here P± = i(l ± 73). To supress the contribution of the region near the 
border of the Brillouin zone we introduce the following interaction of PV 
fields: 

Ipv = I E{^(a:)7M {Pr + KPl) ^(a: + A) - ^x) X 

X [^{x + A) - ^{x)]} + ^ E M^{x)ti^{x) + h.c. (9) 
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Here ^{x) is doublet of bosonic PV fields ip- having the same spinorial 
structure as ijj- and satisfying antiperiodic boundary condition; Pl = 
1(1 + 73T3), Pr = |(1 — 73r3); ri, are Pauli matricies. The regularized 
action 

ICR = ^CW + 

in the formal continuum limit, when the Wilson terms are neglected, is 
invariant under chiral gauge transformations of the fields i/j±, 
The action Icr generates the following propagators: 



iB{p) + W{p) 
B'^{p) + W^{p) 



^ i^(p) + + M l + n iBip)^Wip)-M l-n 

52(j)) + (W(p)+M)'' 2 B^{p)^{W{p)-Mf' 2 ^ 

where ^(p) = E B^{p)j^, B^{p) =B^{p,0) is given by eq.(|5) and is 
given by (|). 

The polarization operator is a sum of the terms which have a form 

N/2 

n,.. - E 9' Tr {^,P±V,{p)S{p)j,P±VMS{p)), 

p=-N/2+l 

where V^{p) is a lattice vertex function and the propagators 5* are given by 
eqs.([TUD, (0). 

Consider firstly the contribution of physical fields. As the model is 
anomaly free, the terms proportional to 73 cansel and one gets 



Following the reasonings of ref . one can show that analogous expression 
for PV fields cancel the contribution of momenta close to the border of the 
Brillouin zone and the remaining part in the limit N ^ 00 coincides with 
the gauge invariant continuum result. Separating the summation domain 
into two parts l^n^ \p\ < 7 < I' Vout- \p\ > and choosing 
M ~ A^^, 5 < 7, one sees that in the domain Kut one can expand the PV 
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fields propagators in terms of M^. Tlie zero order term coincides witli the 
propagator of -i/^-fields and due to our choice of charges the leading terms in 
the expansion of U^^ cancel the value H^J^. The next term is maj orated 
by N^'^^M'^ and vanishes in the limit N ^oo. 

In the domain V^-a one can expand the Wilson term in power of p. The 
leading term in the limit N ^ oo gives the gauge invariant continuum 
expression and the higher order terms are majorated by 7V^t~^ when 

By the same reasoning all higher order diagrams, which correspond to 
power counting convergent integrals in the limit N ^ oo coincide with the 
corresponding continuum expressions. So we expect that simulations of 
the 11112 model (§) regularized by the PV action (0) has to provide a 
reasonable approximation to the continuum toron model. 

A straightforward calculation gives for the regularized lattice 11112 de- 
terminant normalized on 1 at h=0 the following expression: 

DcR = Dcw[h]Dpv[hi (12) 

Dew = Dlw[h]Dl^[2hl 

D+w[h] = n JTT-TTy 

r.i Vt' 2/i]|' + M^[B\p, 2h) + B\p, 0) - 2W\p) + M^lX 

Vpv[ri\ — [[ 

p^-N/2+l 

Here 



H\p, M]H[p,-M] 



G\p, h] = [Boip. h) + i5i(j), h)][Boip. 0) - i5i(j), 0)] + W\p), 

H[p, M] = B\p, 0) + {W{p) + M)\ 

B^{p, h), B^{p, h) and W{p) are defined above. One sees that the regulariza- 
tion changes only modulus of the lattice determinant whereas its argument 
remains intact, i.e. AigDcR = -^TgDcw- 



12 




1,2,3 



■1.0 



-0.5 0.0 



0.5 



1.0 



1.5 h. 



Fig. 4. 11112 determinant arguments AigDc as functions 
of hi at different Hq: 
1 - Arg Dec on the torus; 

2,3- Arg DcR= Arg Dck= Aig Dew' 2 - N = 32, 



3 - 160 



The corresponding expression in the continuum toron 11112 theory looks 
as follows [5-7]: 

Dcc[h] = Dic[h]Dlc[2h]. (13) 



Here 



D+c[h] = e 



iTrhi{ho+ihi) 



Y[F[n,h]F[n,-hi 



n=l 



and expression for the F[n, h] is given above. 

The determinants ([T2|), ([13| ) satisfy the following symmetry properties: 

D[ho, hi] = D*[hi, ho] = D*[-ho, hi] = D*[ho, -hi], 



Dcc[ho, hi] = Dcc[ho + no, hi + ni], no, ni = 0, ±1, ±2, 
DcR[ho, hi] = DcR[ho + A^no, hi + Nm]. 



(14) 
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Fig. 5. 11112 regularized polarization operators IIcr{0)/IIcc{0) 
as functions of M: 
1- N = 160, 2 - 320 

Due to these properties it is sufficient to take the external fields in the 
interval < ho<hi<N/2. 

It follows from ([Tl|) , that only diagrams with more than two external lines 
contribute to the arguments of the lattice and continuum determinants. 
Such diagrams on the lattice differ from the corresponding continuum ones 
by the terms of order 0{1/N). So in the framework of perturbation theory 

ArgDcR = Arg Dew ^ -^^g Dec, when A^^oo. (15) 

Computer simulations of the determinant argument were performed for 
N=32, 160 and are presented at Fig. 4. One sees that eq.([T5]) is valid in 
the interval — 0.5 < /ii < 1.5. 

Since the higher order diagrams are convergent the best agreement for 
the determinant modulus is achieved at the values of M for which the 
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Fig.6. 11112 determinant modulus \Dc\ as functions of hi at ho = 0.2: 

1 - \Dcc\ on the torus; 2, 3, 4 - \Dcr\ with PV field, M = Mq: 

2 - AT = 32, 3 - 160, 4 - 320; 

5, 6 - l^ci^l with counterterm: 5 - = 32, 6 - 160 

regularized polarization operator IIcr{0) is closest to the continuum value 
ncc(0) = 87r. Fig.5 shows that the ratio nci?(0)/ncc(0) for A^>160 differs 
from 1 less than by 2% if M = Mo(iV) - iV"^/l This gives Mo=0.03 - 0.05 
at N=160 - 320. It agrees with the analytical estimates which gives 

Ucr{0) = Stt + 0{1/MN) + OiM^ln^N), if M^O when AT^oo. 

The dependence of the 11112 regularized determinant modulus on hi at 
/io=0.2 for these values of M is given at Fig.6 (curves 1-4). One sees that 
results obtained for the lattice Wilson action with PV regularization agree 
with the continuum theory in a certain interval. When the value grows 
from 32 to 320 this interval expands from |/ii|<0.25 to — 0.5</ii<0.9. 

As in the case of vectorial theory another way to achieve the agreement 
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of lattice and continuum 11112 theories is to introduce to the action (@) 
a real gauge noninvariant counterterm. Then we get for the 11112 lattice 
determinant 

DcK = Dew exp [kcihl + hj)] , 

where Dew defined above. Choosing /c(7=6.4005 we get the curves 5, 6 
shown at Fig. 6. At N=160 one has an agreement of the theory with coun- 
terterm and continuum one in the interval — 0.5 < /ii < 1.5. 

4 Discussion 

In this paper we showed that introducing to the fermion lattice action 
additional Pauli - Villars type regularization one can supress the gauge 
symmetry breaking effects caused by the Wilson term both in perturba- 
tive and nonperturbative regime. No gauge noninvariant counterterms are 
needed to get the correct continuum result. For a finite lattice spacing there 
are symmetry breaking effects of order a. In principle one can avoid the 
symmetry breaking completely by introducing a manifestly chiral gauge in- 
variant action like SLAC model supplemented by PV regularization. Our 
calculations show that in the vectorial model it produces a manifestly gauge 
invariant result in a good agreement with the continuum expression. How- 
ever in the anomaly free chiral models there are some problems related to 
the nonlocality of the SLAC model. These problems will be discussed in a 
separate publication. 
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